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摘 要 : 在 塔里木 河 下 游 沙漠 地 区 防护 林 建 设 过 程 中 ,地 下 成 水 灌溉 条 件 下 土壤 剖 
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显 , 且 存在 一 定 的 规律 性 。 为 掌握 土壤 垂 向 水 盐 动 态 ,建立 数值 模型 是 十 分 高 效 的 手段 ,通过 重建 
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续 可 持续 灌溉 管理 。 通 过 HYDRUS-1D 模 拟 防 护林 建设 初期 不 同 矿 化 度 地 下 水 灌溉 条 件 下 土壤 水 
盐分 布 特征 和 时 间 变 化 ,并 根据 实测 的 土壤 水 盐 数 据 , 优 化 土壤 水 分 特征 参数 ,评价 模型 的 适用 
性 。 结 果 表明 :表层 (0~30 cm) 土 壤 含 水 量 和 含 盐 量 受 灌溉 影响 大 ,数值 波动 剧烈 , 且 土 壤 盐 分 表 


聚 强烈 ;深层 (50~150 cm) 则 受 灌溉 的 影响 小 ,数值 波动 小 。 数 值 模型 模拟 值 与 实测 值 吻 合 ,能 较 
好 地 反映 土壤 水 盐 运 移 情况 ,建立 的 模型 可 用 于 模拟 不 同 矿 化 度 灌 溉 条 件 下 塔里木 河 下 游 流动 沙 
漠 地 区 土壤 水 盐 运 移 状况 ,为 该 地 区 生态 防护 林 建 设 提供 理论 基础 和 技术 支撑 。 

X 键 词 : 成 水 灌 溉 ; 土壤 水 盐 运 移 ; HYDRUS-1D 模型 ; 塔里木 河 下 游 
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格 ( 尔 木 ) 库 ( 尔 勒 ) 铁 路 是 中 国 通 往 西亚 `. 地 中 
海 和 黑海 地 区 的 陆运 大 通道 ,对 国民 经 济 发 展 和 巩 
回国 防 意义 重大 。 铁 路 穿越 塔里木 河 下 游 生态 脆 
弱 区 ,风沙 活动 强烈 ,需要 建设 生态 防护 林 以 确 
保 铁 路 安全 运行 ,但 由 于 该 区 域 水 资源 尤其 是 地 表 
水 优先 用 于 农业 生产 。 而 地 下 水 受 蒸发 和 排泄 过 
程 影 响 , 矿 化 度 很 高 ,几乎 无 法 直接 用 于 常规 的 植 
物 灌溉 ,利用 这 种 劣质 的 地 下 水 灌溉 超 耐 盐 灌 木 ， 
是 该 地 区 防护 林 建 设 的 重要 技术 途径 。 研 究 表明 ， 
咸 水 灌 溉 对 大 多 数 荡 漠 植 物 的 生长 有 显著 的 抑制 
作用 ,高 盐水 灌溉 更 会 导致 植物 死亡 ”。 目 前 ， 
内 外 学 者 已 经 针对 咸 水 灌 溉 做 了 大 量 研究 ,内 容 
涉及 土壤 物理 .化 学 和 生物 ,也 包含 对 土壤 水 盐 
运 移 及 盐分 分 布 的 分 析 ” 。 但 以 往 研究 主要 集中 
在 矿 化 度 小 于 5 gL 的 咸 水 , 最 大 不 超过 30 gL, 
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对 于 30 gL' 以 上 咸 水 灌 溉 的 研究 ,几乎 没有 涉及 ， 
而 高 盐 灌溉 环境 下 土壤 的 粘 清 系 数 .弥散 度 均 在 一 
定 程度 上 增加 ,土壤 水 盐 运 移 规 律 与 低 盐 情况 也 会 
有 不 同 。 而 根据 实测 , 拟 建 防护 林 区 地 下 水 TDS( 溶 
解 性 固体 总 量 ) 普 遍 超过 30 gL', 最 高 时 可 达 50~ 
70 gL 以 上 。 因 此 ,从 地 下 水 利用 的 角度 考虑 , 探 
讨 土 壤 放 面 水 盐分 布 特征 ,对 荒漠 区 地 下 成 水 资源 
的 高 效 利 用 及 耐 盐 碱 车 演 植物 的 栽培 管理 ,具有 深 
建立 水 盐 运 移 模 型 是 分 析 土 壤 水 盐 运 移 的 重 
要 手段 之 一 , 相 比 于 传统 的 野外 测量 ,使 用 模型 模 
拟 水 盐 规 律 能 够 节省 大 量 的 人 力 物 力 。 目 前 ,国内 
外 广泛 应 用 的 几 种 土壤 水 盐 运 移 模型 有 HYDRUS、 
PROFLOW WAV-ES .2DFATMIC SHAW 等 。 其 中 ， 
HYDRUS-1D 模型 是 由 美国 盐 改 中 心 开 发 的 免费 用 
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于 模拟 变 饱 和 多 孔 介 质 水 分 .热量 溶质 迁移 的 模 
型 程序 ” ,能 够 结合 大 气 过 程 ( 蒸 散 、 降 十 灌溉 )、 土 
壤 水 分 及 溶质 运 移 .热量 .植物 根系 吸水 等 过 程 进 
行 综合 模拟 。 近 年 来 ,HYDRUS 软件 被 广泛 应 用 于 
模拟 不 同 土地 利用 方式 下 的 土壤 水 分 运动 ”中 和 溶 
质 运 移 "等 ,并 被 证 明 有 很 好 的 适用 性 和 有 效 性 。 

本 文 以 不 同 矿 化 度 地 下 水 滴 汉 条件 下 ,夏秋 两 
季 各 一 个 完整 灌水 周期 的 土壤 水 分 实测 数据 以 及 
夏秋 两 季 灌 溉 周期 末 土 壤 齐 面 盐 分 数据 为 基础 , 结 
合 当地 的 气象 .灌溉 、 植 被 情况 ,运用 HYDRUS-ID 
FRAI XT RE (Tamarix chinensis Lour. ) 滴 灌 带 内 土壤 
水 盐 运 移 过 程 进行 了 模拟 ,通过 对 模拟 值 与 实测 值 
比较 分 析 , 验 证 了 模型 的 可 靠 性 和 精度 。 人 研究 结果 
可 以 为 该 区 域 防 护林 工程 建设 提供 理论 依据 和 有 
效 的 技术 文 撑 。 


1 数据 与 方法 


1.1 研究 区 概况 

研究 区 位 于 塔里木 盆地 东北 缘 农 二 师 三 十 三 
团 格 库 铁 路 沿线 沙漠 中 ,地理 位 置 40"45'N , 87°18'E 
(图 1)。 该 地 区 年 平均 气温 12 % ,极端 最 高 气温 可 
达 44.3 % ,极端 最 低 气 温 则 超过 -30.9 C ,年 平均 降 
水 量 24.7~54.1 mm, 年 潜在 蒸发 量 2000~3000 mm, 
每 年 4 一 7 月 为 风 季 ,主导 风向 以 东风 ,东北 风 为 主 ,年 
平均 风速 1.6~3.2 mes ,最 大 风速 13.3~20.5 mest, ik 
验 区 土壤 质地 属于 砂 士 ,试验 期 间 地 下 水 埋 深 在 
4.25~4.65 m 之 间 波 动 。 
1.2 试验 方法 

本 文采 用 田间 控制 试验 ,在 树 柳 种 植 试验 区 林 
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1 研究 区 示意 图 
Fig. 1 Geographical location of the study area 
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地 内 进行 。 试 验 区 设 有 水 源 井 ,用 于 地 下 水 位 监测 
与 抽水 灌溉 。 为 研究 不 同 高 矿 化 度 咸 水 灌 溉 对 植 
物 和 土壤 的 影响 ,对 当地 地 下 水 进行 淡化 处 理 , 使 
用 滴灌 造林 。 灌 水 矿 化 度 梯度 分 4 个 (10 g L',25 
g L',40 g: ,原水 ) ,其 中 原水 指 无 处 理 地 下 水 。 
灌溉 水 的 矿 化 度 稳定 性 如 图 2 所 示 , 除 原水 变化 幅 
度 较 大 外 (33.32~54.14 g L) ,淡化 设备 基本 可 以 将 
灌溉 水 稳定 性 维持 在 图 2 所 示 的 矿 化 度 水 平 上 。 
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图 2 试验 期 间 灌溉 水 矿 化 度 监 测 
Fig.2 Monitoring of irrigation water salinity during 


the experiment 


野外 监测 时 间 为 2019 年 5 月 30 日 至 2019 年 9 
月 20 日 。 监 测 指标 包括 树 柳 根系 分 布 、 株 高 . 叶 面 
积 指数 .土壤 含水 量 .土壤 电导 率 (EC,。) 地 下 水 位 、 
地 下 水 电导 率 和 气象 数据 。 

根系 密度 的 空间 分 布 采 用 生长 季 末 根系 整体 
挖掘 取 土 , 洗 根 烘 干 分 层 称 重 来 确定 ,重复 3 次 。 采 
用 中 子 土 壤 水 分 仪 测定 土壤 水 分 含量 ,将 中 子 仪 测 
点 置 于 滴 头 下 方 ,监测 深度 范围 为 0~100 cm( 本 文 
选用 土 层 深度 分 别 为 15 cm 30 cm 50 cm .90 cm 处 
的 土壤 水 分 含量 作为 本 文 模型 验证 数据 )。 土 壤 电 
导 率 通过 土 外 分 层 取样 (测定 深度 为 0~160 cm,20 
cm 一 层 ),1:5 土 水 比 测定 (ECi) ,使 用 Excel 拟 合 公 
式 将 电导 率 值 转化 为 含 盐 量 ,8 月 1 日 和 9 月 20 日 (2 
次 灌溉 周期 末 ) 为 取样 时 间 ,对 土壤 进行 分 层 取 样 ， 
重复 3 次 。 滴 灌 带 滴 头 间距 为 60 cm, 滴 头 流量 为 
3.8 Leh! ,滴灌 带 间 距 2.5 m , 滴灌 带 长 度 可 达 100 m, 
单 次 连续 灌水 8h ,滴灌 湿润 宽度 可 达 65 cmo 
1.3 数值 模拟 

使 用 HYDRUS-1D 软件 模拟 2019 年 的 土壤 水 盐 
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运 移 情况 ,模拟 土 层 深度 范围 为 0~160 cm。 模 拟 时 
段 从 2019 年 5 月 30 日 至 9 月 30 日 ,共计 124 d, 使 用 
7 月 22 日 至 8 月 1 日 和 9 月 10 日 至 9 月 20 日 2 个 灌 
水 周期 的 实测 土壤 水 分 数据 .8 月 1 日 和 9 月 20 日 士 
壤 剖 面 盐分 含量 实测 数据 对 模型 进行 验证 。 初 始 
时 间 步 长 设 为 0.001 d, 最 大 时 间 步 长 为 0.5 do RK 
迭代 次 数 设 置 为 10, 允许 程序 在 时 间 步 长 范围 内 对 
时 间 步 长 自动 调整 。 

采用 Van Genuchten-Mualem 模型 作为 土壤 水 分 
特征 曲线 模型 。 将 治水 周期 前 测定 的 土壤 水 分 、 电 
导 率 作为 初始 条 件 输入 模型 。 根 据 当 地 的 气候 情 
况 灌水 情况 ,确定 模拟 土 体 的 上 边界 为 开放 的 大 
气 边界 ,接受 灌溉 补给 、 土 面 区 发 .降水 等 ,潜在 落 
发 量 ET, 根 据 FAO 推 荐 的 Penman-Monteith 公 式 由 气 
象 数据 和 叶 面 积 指数 计算 得 到 。 研 究 区 地 下 水 位 在 
4.25~4.65 m 之 间 浮 动 , 根 据 周 宏 ASCE AF 
究 , 砂 土 的 毛管 水 上 升 高 度 在 潜水 层 以 上 120 em, 
因而 地 下 水 对 模拟 土 体 (0~160 em) 几乎 无 影响 , 故 
下 边界 条 件 设 定 为 自由 排水 。 
1.3.1 水 分 运动 模型 ”本 次 模拟 运用 Richards 方 程 
描述 变 饱和 区 一 维 垂直 方向 上 的 土壤 水 分 运动 : 


ils AKOE + n| ~ S(z, 0) (1) 


式 中 :6 为 体积 含水 量 (cm +em™) ;为 时 间 (d);z 是 垂 

直方 向 的 土壤 深度 坐标 (cm) , 取 地 表 为 原点 ,向 下 为 

TEs K(h) AK 8 (em d”) sh Ay RIK (om) 5S 为 源 

汇 项 ,主要 用 来 描述 根系 吸水 情况 (cm sem? +d"). 
非 饱和 土壤 水 力学 参数 描述 如 下 : 


0 = Ue h<0 
6( 和 =1 — (1 +|ah|') (2) 
6,,h>0 
K(h)=K Sji- (1-3) (3) 


SUH 00,2 EEF HL AES AK et RAS EKE Cem- 
cm ) ;a mn 为 经 验 参 数 , 其 中 m=1-1/n;Ks 为 饱和 
导 水 率 (cm*d"');5. 为 无 量 纲 的 相对 含水 量 ;/ 为 和 孔 际 
连接 系数 ,通常 取 0.5。 

1.3.2 盐分 运动 模型 ”在 不 考虑 固 相 吸附 的 情况 
下 ,选用 平衡 态 的 对 流 弥 散 方 程 来 描述 根 区 的 溶质 


运 移 : 


aC _ a oC\_ v6C 
ot = Ż(oD c) dz (4) 


SOA: CEN IK AYRE (mo. cem’); D HANK BAK 


PAB Cem d”) o RIR PIILEE R cmd”) 
1.3.3 根系 吸水 模型 ， 因 潜在 蒸腾 量 ( 即 潜在 根系 
吸水 量 ) 在 根 区 的 分 布 情况 与 根系 密度 的 空间 分 布 
相 一 致 ,本 文 根 据 实地 分 层 挖掘 根系 的 结果 ,将 各 
层 根 系 干 重 输入 模型 ,由 模型 转化 生成 根系 分 布 函 
数 。 模 型 将 根系 吸水 参数 转化 为 水 分 胁迫 函数 来 
表示 根系 吸水 情况 ,本 文采 用 Feddes 梯形 函数 来 描 
述 树 柳 根系 吸水 过 程 ,根据 Moayyad 和 Crinevskii 的 
研究 结果 ,其 吸水 参数 取 值 见 表 1。HYDRUS-1D 
中 的 根系 吸水 采用 修正 的 Feddes 模 型 ,如 式 5 所 示 : 

S(h, h,, x, z)=a(h, h,, x, z)b(x, z)S,T, (5) 
式 中 : alh, h, x, 2) H ERKENE PAK; h, 为 渗透 
JEH (cm); bl, z) 为 根系 分 布 函数 (cem”) ;$, 为 与 蒸腾 相 
关 的 地 表 长 度 (cm) ;7 为 潜在 腾 发 速率 (cmd')。 


表 1 根系 吸水 参数 


Tab. 1 Water absorption parameters of roots 


hiem hem jcem hem hyem Paem':d | ra/em*d" 
-0.1 -2.0 -80.0 -250.0 -15000.0 0.5 0.1 
注 : 记 为 根系 吸水 厌 氧 点 土壤 基质 势 ; 户 为 根系 吸水 最 适 点 开始 土 
壤 基 质 势 ; 户 为 根系 吸水 最 适 点 土壤 基质 势 (ma 为 -80 土壤 基质 势 ， 
有 ;为 -250 土 壤 基 质 势 ) ;为 根系 吸水 萎 薄 点 土壤 基质 势 ;ry 和 5 
分 别 为 2 个 假设 的 潜在 作物 蒸腾 率 。 


1.3.4 潜在 蕉 散 量 的 计算 ”气象 数据 作为 上 边界 条 
件 输入 模型 时 ,采用 Penman-Monteith 公式 计算 潜在 
REE ,该 公式 在 干旱 环境 下 的 适用 性 已 得 到 许多 
学 者 的 验证 2 ,其 主要 公式 为 : 
_0.4084(R, -6) + ya oe) 
° A+y(1+0.34u,) 

RE : ET HRPE (mm: h’); R, WEAR A EIT 
[LMJ md”];C 为 土壤 热 通 量 LMJ (md) ]; yF 
湿 表 常数 ;7T 为 2 m 高 处 日 平均 气温 (%C);ww 为 2m 高 
处 日 平均 风速 (ms ie. 为 饱和 水 气压 (kPa) ;e, 为 
实际 水 汽 压 (kPa) ;A 为 饱和 水 汽 压 随 温度 的 斜率 。 
1.4 模型 的 校准 与 评价 

根据 实测 土壤 质地 的 垂直 分 布 特征 ,将 不 同 深 
度 的 土壤 谢 面 分 为 5 层 (具体 土 层 深度 划分 见 表 
2)。 各 层 水 分 特征 参数 由 HYDRUS-1D 模 型 的 Ros- 
seta 模 块 , 输 入 各 层 实 测 土壤 粒 径 和 容重 ,得 到 各 层 
的 初始 参数 ,使 用 夏季 灌水 周期 内 实测 土壤 水 分 数 
据 进 行 参 数 校准 ,土壤 纵向 弥散 度 为 60 cm。 为 提 
高 数值 模拟 的 准确 性 ,需要 对 土壤 水 分 特征 参数 进 


(6) 
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表 2 研究 区 土壤 粒 径 及 土壤 水 分 特征 参数 


Tab.2 Soil particle size and soil water characteristic parameters in the study area 


土 层 颗粒 组 成 土壤 容重 土壤 水 分 特征 参数 
深度 /cm 砂粒 粉 粒 。 BAL /gcm bcm cm’ OQ./cm tem” 天 /cmd alem” n l 
0~20 79.6 19.6 0.8 1.46 0.015 0.382 100 0.047 2.82 0.5 
20~60 98.9 1.1 0.0 1.52 0.013 0.377 56 0.031 2.59 0.5 
60~80 97.5 25 0.0 1:55 0.011 0.368 86 0.032 3.26 0.5 
80~100 94.1 5.9 0.0 1.58 0.017 0.370 56 0.035 3.13 0.5 
100~150 97.5 2:5 0.0 1.54 0.013 0.372 100 0.032 4.04 0.5 


行 调整 。 根 据 敏 感性 分 析 , 主 要 对 饱和 导 水 率 Ks 和 
经 验 参数 nn 进行 优化 ,优化 后 的 结果 见 表 2。 将 水 盐 
观测 数据 分 为 两 段 ,2019 年 7 月 22 日 到 2019 年 8 月 
1 日 用 于 模型 校准 , 另 一 时 段 用 于 模型 验证 。 采 用 
均 方 根 误差 (RMSE) .决定 系数 ( 尼 ) 和 模型 效率 系数 
(NSE )3 种 方法 来 评价 本 文 数值 模拟 的 效果 。 


2 结果 与 分 析 


2.1 土壤 水 分 模拟 
从 图 3 可 以 看 出 ,土壤 水 分 的 模拟 值 与 实测 值 


0.25- (a) 10 g LK 


一 ”15 cm 模拟 值 e15 cm 实测 值 
z i ----30cm 模 拟 值 ”o 30 cm 实测 值 
B EMI BN ge 50 cm 模拟 值 a 50 cm 实测 值 
加 ~…… 90 cm 模拟 值 A 4 90 cm 实测 值 
© 0.15 
ist 
A 0.10 
x 
Hoste) NA 

OOD a a a la aa a 
cy D N SN D 9 
SPP LP KL LK S 
日 期 /月 -日 
(c) 40 gL 灌溉 水 

» 0.20 一 一 15 cm 模拟 值 。 。15 cm 实测 值 
H o ---- 30 cm 模拟 值  o30cm 实 测 值 
5 ae 50 cm 模拟 值 。 450 cm 实测 值 
g 0.15 从 sree 90 cm 模拟 值 。 ^90 cm 实测 值 
At 
x 
a 
ES 
H 


E ORAS EKE 5 0AE IKE KARTKE 5 EL 3 n HS FLBR ACN Sd RARR AIRERA 


比较 吻合 。 灌 溉 后 ,表层 (0~30 em) 土壤 含水 率先 
迅速 增加 ,达到 15% ~20% ,然后 快速 下 降 10% 以 
下 。50~90 cm 土 层 含水 量变 化 较 小 , 土 层 含水 量 通 
党 在 5%~10% 之 间 波 动 , 其 中 土 层 90 cm 处 对 灌溉 
的 响应 有 比较 明显 滞后 ,数值 变化 不 是 很 剧烈 。 水 
分 再 分 布 阶段 ,表层 十 壤 含水 率 低 于 下 层 , 说 明 表 
层 土壤 受 灌溉 入 渗 和 莱 散 的 影响 剧烈 ,模型 能 够 较 
好 地 模拟 出 土壤 水 分 从 入 渗 到 再 分 布 过 程 的 变化 
情况 。 以 上 结果 说 明 ,HYDRUS-1D 模 型 对 研究 区 
土壤 水 分 运 移 的 模拟 适用 性 较 强 。 


0.20r (b) 25 gL Wk 


7 一 一 15 cm 模拟 值 e 15 cm 实测 值 
h ---- 30 cm 模拟 值 。 0 30 cm 实测 值 
2 0.15 + Bee 50 cm 模拟 值 '' 450 cm 实测 值 
g Iv 一 90 cm 模拟 值 ，\，^90 cm 实测 值 
® ool pha N 
x 
如 
H 0.05 上 È 
fe} 
0.00 1 5 1 L 74 L 1 © J 
py P o S OOA 
TEE REF: 
日 期 /月 -日 
d) 原水 
0.15 r @ 一 一 15 cm 模拟 值 。 。15 cm 实测 值 
z o ---- 30 cm 模拟 值  o30 cm 实测 值 
9 a 二 50 cm 模拟 值 450 cm 实测 值 
$ Ao. core 90 cm 模拟 值 。 ^90 cm 实测 值 
w 
x 
<a 
EN 
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图 3 土壤 含水 率 模拟 值 与 实测 值 


Fig.3 Simulated and measured values of soil moisture content 
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对 土壤 含水 率 总 体 模拟 效果 进行 评价 ,人 研究 结 
果 表 明 经 过 校准 的 HYDRUS- OR Aa AET 
的 模拟 土壤 水 分 动态 变化 (图 3) ,这 主要 体现 在 较 
小 的 RMSE 、 较 高 的 尾 和 NSE。 通常 党 而 言 ,RMSE 表 
示 模 拟 值 与 实测 值 之 间 的 偏差 ,但 不 能 反映 模拟 与 


培训 面 盐 分 模拟 值 与 实测 值 的 对 比 ,模拟 值 与 实测 
值 基 本 吻合 。 从 图 4 可 以 看 出 ,在 不 同 灌溉 水 矿 化 
度 处 理 下 ,土壤 含 盐 量 从 表层 到 深层 呈现 先 快速 降 
低 ,后 保持 稳定 的 分 布 特征 。 土 壤 表层 存在 较为 明 
显 的 积 盐 区 , 随 着 灌水 矿 化 度 的 增加 ,各 层 土壤 含 


实测 的 拟 合 情况 ,而 许 和 NSE 却 能 很 好 地 解释 模型 

以 合 效 果 。4 种 矿 化 度 灌溉 下 土壤 水 分 RMSE 在 
0.014~0.020 之 间 , 尼 在 0.71~0.88 之 间 ,NSE 在 0.67~ 
0.85 之 间 。 

在 设置 的 4 种 矿 化 度 灌溉 下 ,不 同 土 层 的 评价 
效果 见 表 3。 不 同 深度 土壤 水 分 在 校准 和 验证 期 间 
AY RMSE 值 均 在 0.005~0.032 之 间 ,从 数值 角度 看 ， 
模型 模拟 效果 较为 理想 ,具有 一 定 的 可 靠 性 ;而 尼 
和 NSE 这 2 种 评价 之 间 反 馈 的 结果 ,不 同 土 层 模拟 
效果 为 :30 cm>15 cm>50 cm>90 cm,50 cm 和 90 cm 
土 层 模拟 精度 低 ,甚至 在 90 cm 土 层 的 模型 效率 系 
数 显 示 为 不 可 信 (NSE<0)。 表 层 土壤 (0~30 cm) 的 
模拟 效果 要 好 于 下 层 (50~90 cm). 

2.2 土壤 盐分 模拟 
图 4 为 2019 年 8 月 1 日 和 2019 年 9 月 20 日 的 土 


盐 量 也 相应 增加 ,土壤 表 聚 也 越 明 显 。 

由 表 4 可 知 , 在 4 种 矿 化 度 灌溉 下 ,土壤 盐分 模 
拟 值 与 实测 值 RMSE 分 别 为 0.13 ,0.23 .0.33 0.59, R? 
分 别 为 0.92 .0.85 , 0.86 Fil 0.74, NSE 分 别 为 0.92、 
0.78、0.89 Fl 0.65. #&/)\A9 RMSE, #8 fey AY R Al NSE 
表明 不 同 矿 化 度 灌 溉 下 ,土壤 盐分 模拟 值 与 实测 值 
的 空间 分 布 规律 较为 一 致 ,差异 性 不 大 。 
2.3 情景 模拟 

为 了 更 好 地 反应 未 来 时 段 咸 水 灌溉 后 ,林地 土 
壤 水 盐 运 移 过 程 的 变化 趋势 ,本 文 设 定 了 下 一 年 度 
灌溉 起 点 至 植物 生长 季 末 (5 月 30 日 一 9 月 22 日 ) 约 
4 个 月 的 模拟 时 长 。 

水 分 特征 参数 采用 调 参 后 的 夏季 灌溉 参数 , 选 
中 HYDRUS-1D 中 植物 生长 模块 ,将 生长 期 间 监 测 记 
录 的 树 柳 株 高 . 叶 面 积 指数 、 根 系 深 度数 据 输 入 模 


表 3 土壤 含水 率 模 拟 结果 统计 分 析 


Tab.3 RMSE, R’ and NSE of simulated value and measured value of soil moisture content 


灌水 | ca RMSE | R E NSE E 
矿 化 度 /g*L" 校准 验证 校准 验证 校准 验证 
10 15 0.027 0.021 0.83 0.83 0.75 0.79 
30 0.020 0.013 0.83 0.91 0.50 0.85 
50 0.026 0.021 0.54 0.58 0.29 0.07 
90 0.025 0.012 0.36 0.58 -5.20 -0.04 

总 体 效 果 0.020 0.75 0.67 
25 15 0.021 0.026 0.79 0.69 0.66 0.72 
30 0.023 0.030 0.81 0.72 0.75 0.64 
50 0.024 0.023 0.83 0.88 0.75 0.83 
90 0.017 0.016 0.62 0.79 0.47 0.71 

总 体 效 果 0.020 0.77 0.72 
40 15 0.015 0.011 0.91 0.91 0.82 0.81 
30 0.021 0.008 0.82 0.93 0.59 0.91 
50 0.032 0.013 0.81 0.71 0.67 0.79 
90 0.019 0.012 0.58 0.53 -0.60 -0.13 

总 体 效果 0.014 0.88 0.85 
原水 15 0.022 0.007 0.75 0.93 0.68 0.89 
30 0.012 0.007 0.83 0.85 0.77 0.76 
50 0.022 0.009 0.69 0.67 0.63 0.69 
90 0.019 0.005 0.40 0.52 -3.48 -3.82 

总 体 效果 0.015 0.71 0.68 


TE: RMSE 为 均 方 根 误差 ; 尼 为 决定 系数 ;NSE 为 模型 效率 系数 。 下 同 。 
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(al) 10 gL“ YE 7K (2019-08-01) (a2) 10 gL -灌溉 水 (2019-09-20) 


土壤 含 盐 量 /g.kg-1 土壤 含 盐 量 /g.kg-! 
0.5 10 1.5 20 2.5 


g 0 0.51.0 15 2.0 2.5 g 0 

= 40 " = 40 : 

me 30 me 80 。 。 实测 值 
+4 120 H 120 一 模拟 值 


(b1) 25 g:L- 灌 溉 水 (2019_-08-00) (b2) 25 g. 志 灌溉 水 (2019_09_20) 


土壤 含 盐 量 /g.kgr7 土壤 含 盐 量 /g.kg-! 
41S 20 25 30 015 20 2.5 3.0 3.5 
g T T 3 1 fz] T T T Pal 
R OF ig 40 
ie - 实测 值 E 2 。 SE 
H 120 一 模拟 值 H 120[ | ”一 模拟 什 
160 160L " 
(cl) 40 g-L“#E387kK (2019-08-01) (c2) 40 gL-! 浇 溉 水 (2019-09-20) 
土壤 舍 盐 量 /g-kg! 土壤 含 盐 量 /g.kgr 
4 0 3 4 5 7 a a 5 6 
2 80 2 
E all © 实测 值 ak 
H " 一 模拟 值 4 
1604" 
(d1) 原水 (2019-08-01) (d2) 原水 (2019-09-20) 
土壤 含 盐 量 /g.kgr 土壤 含 盐 量 /g.Kg-+ 
0 3 4 5 6 7 0 3 4 5 6 7 
9 T T T T 1 9 T T T T 1 
R 40 . R 40 
KE 80| 叫 " u 实测 值 KK 80; «lf, a 实测 值 
当 120 ， 一 模拟 值 a ol} ”一 模拟 什 
160 160 


图 4 土壤 盐分 模拟 值 与 实测 值 对 比 
Fig. 4 Comparison of simulated and measured values 


of soil salinity 


表 4 土壤 盐分 模拟 值 与 实测 值 的 RMSE, R* 和 NSE 
Tab. 4 RMSE, R’ and NSE of simulated and 


measured soil salinity values 


ee ee es 
RMSE 0.13 0.23 0.33 0.59 

R 0.92 0.85 0.86 0.74 
NSE 0.92 0.78 0.89 0.65 


型 ,并 默认 从 灌溉 起 点 到 生长 季 末 期 植物 的 根系 生 
长 是 均匀 的 。 结 合 实际 灌溉 数据 和 逐日 气象 数据 ， 
对 土壤 水 盐 进 行 模拟 ; 由 于 试验 地 植物 的 根系 主要 
分 布 在 0~40 cm 范围 , 且 上 文中 土壤 水 分 在 此 区 间 
模拟 精度 较 高 ,因此 主要 模拟 15 cm 和 30 cm 深度 土 
JZ ,模拟 结果 如 图 5 所 示 。 

在 经 过 灌溉 后 ,土壤 盐分 和 水 分 是 有 紧密 联系 
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的 ,尤其 以 表层 土壤 最 为 明显 。 如 图 5a 和 图 5b 表 示 
土壤 含水 率 变 化 过 程 ,其 中 每 次 含水 量 陡 升 表示 在 
该 天 进行 了 灌溉 (降水 较 少 ,可 忽略 )。 如 图 5c 和 图 
5d, 随 着 灌溉 开始 后 ,土壤 中 的 盐分 溶解 于 水 中 ,而 
水 分 受 重力 作用 向 下 迁移 ,盐分 随 着 水 分 淋 溶 到 深 
层 土壤 中 ,这 一 过 程 称 为 洗 盐 过 程 , 此 时 土壤 盐分 
快速 降低 。 济 溉 过 程 结束 后 ,由 于 荒漠 地 区 燕 发 强 
烈 ,因此 表层 水 快速 蒸发 ,下 层 土壤 水 分 携带 着 盐 
分 向 表层 迁移 聚集 ,造成 表层 盐分 积累 。 且 随时 间 
推移 ,土壤 水 分 经 过 向 上 运动 和 向 下 运动 ,水 分 含 
量 减 少 , 表 层 土 壤 的 盐分 增加 速率 也 逐渐 变 缓 。 在 
分 析 表 层 土壤 含水 率 和 含 盐 量 在 某 一 周期 内 变化 
过 程 ,二 者 之 间 动 态 变化 关系 密切 。 而 随 着 试验 区 
进行 多 次 灌溉 后 ,除了 原水 灌溉 由 于 矿 化 度 变 化 较 
大 ,其 他 矿 化 度 治 溉 下 表层 土壤 出 现 盐 分 累积 增加 
现象 。 


3 讨论 


塔里木 河 下 游 沙 漠 区 气候 干旱 少雨 ,自然 环境 
恶劣 ,应 用 沙 地 戌 水 滴灌 技术 建设 防护 林 , 能 够 为 
格 库 铁 路 安全 运行 提供 强 有 力 的 保障 。 本 文通 过 
10 g*L'、25 g- L'.40 gL"! 和 原水 (无 处 理 地 下 水 )4 
种 矿 化 度 灌 溉 水 对 土壤 水 盐 运 移 影响 进行 了 人 研究， 
发 现在 灌水 结束 后 ,土壤 表层 含水 率 会 迅速 地 降低 
到 10% 以 下 ,这 与 丁 新 原 等 ” 提 到 的 现象 相 一 致 ; 
水 分 是 保障 植物 正常 生长 的 关键 因素 ,防护林 沙 质 
土壤 的 某 水 性 较 差 ,如何 延 缓 灌水 后 浅 层 土壤 水 分 
的 迅速 流失 ,是 以 后 工作 的 重点 。 

土壤 盐 渍 化 是 防护 林 后 期 建设 与 维护 过 程 中 
经 常 遇 到 的 问题 。 本 文通 过 夏季 和 秋季 2 次 盐分 观 
测 的 结果 ,发 现 灌溉 水 矿 化 度 越 高 ,土壤 表 聚 现象 
也 就 越 明显 ,这 也 侧面 印证 了 窦 山 歌 等 的 结论 ; 
周智 彬 等 六 通过 观测 4 个 年 份 的 水 盐 运 动 规律 , 认 
为 在 灌水 和 营 发 主导 下 ,土壤 整体 含 盐 量 可 在 短期 
内 维持 一 个 相对 平衡 状态 ,并 呈现 出 周期 性 的 波动 
变化 ,本 文 运 用 模型 模拟 还 原 整个 植物 生长 季 根 系 
层 土壤 水 盐 波动 情况 ,也 显示 了 土壤 盐分 具有 周期 
性 的 波动 变化 特点 。 然 而 ,从 长 远 角 度 看 ,在 连续 
进行 高 矿 化 度 灌水 的 影响 下 ,土壤 盐分 源源 不 断 地 
输入 土壤 ,在 各 个 深度 均 有 绥 慢 累积 的 趋势 , 随 着 
苗木 逐渐 成 长 ,根系 的 聚 盐 效应 3 也 开始 显现 ,如 
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图 5 土壤 15 cm 30 cm 剖面 水 盐 运 移 过 程 模 拟 


Fig. 5 Simulation of soil water and salt migration at depths of 15 cm and 30 cm 


何平 衡 和 调控 盐分 空间 分 布 ,并 寻找 有 效 的 脱盐 途 
径 ,是 下 一 步 研究 的 方向 。 

近年 来 ,运用 软件 模拟 水 盐 动态 的 方式 来 参 
与 水 盐 管 理 调控 工作 ,成 为 逐渐 兴起 的 方式 。 余 
He WSS) 、 李 亮 等 乓 使 用 HYDRUS 模 型 对 内 蒙古 河套 
干旱 区 水 盐 运 移 规律 进行 模拟 研究 ,为 区 域 节 水 控 
盐 和 灌区 长 效 发 展 提供 了 新 思路 ; 李 开 明 "5 、 石 培 
君 等 ”运用 HYDRUS 模 型 模拟 了 明 沟 、 暗 管 排水 条 
件 下 土壤 水 盐 运 移 规 律 ,证 明 在 此 条 件 下 ,HYDRUS 
模型 能 够 有 效 地 得 到 应 用 ; 李 会 杰 ”\ 王 在 敏 等 3 
使 用 HYDRUS 模型 分 别 模拟 了 玉米 根 区 水 分 规律 
和 棉花 生育 期 水 盐 运 移 规 律 ,表明 模型 能 够 用 于 在 
作物 灌溉 过 程 中 对 土壤 水 盐 状况 的 监测 。 根 据 以 
上 研究 结果 ,HYDRUS 模 型 能 够 在 多 种 条 件 下 灵活 
应 用 ,并 被 证 明 具 有 相当 的 可 靠 性 和 稳定 性 。 运 用 
HYDRUS 数 值 模拟 手段 ,可 以 在 有 效 降 低 人 研究 成 本 
的 情境 下 ,实现 对 未 来 土壤 水 盐 环 境 变 化 较为 准确 
的 预 佑 ,对 于 其 他 土壤 过 程 的 研究 也 能 起 到 有 效 的 
辅助 作用 。 

在 高 矿 化 度 咸 水 持 续 灌 溉 的 严酷 条 件 下 ,尤其 
是 灌溉 水 矿 化 度 超过 30 g:L7' 时 ,土壤 中 的 水 体 在 


筷 际 介质 中 的 流动 情况 受 土 壤 吸 附 力 、 土 体 孔 际 结 
构 改 变 等 情况 的 影响 ,土壤 溶质 运 移 参数 会 随 之 发 
生 改 变 ;植物 在 高 矿 化 度 咸 水 灌溉 中 受 水 分 和 盐分 
双重 胁迫 的 影响 ,生理 生态 方面 发 生 的 显著 变化 ， 
这 些 问题 都 或 待 进一步 开展 研究 工作 。 滴 灌 后 水 
分 的 实际 运动 是 三 维 的 ,而 本 文 主要 考虑 的 是 滴 头 
下 方 水 分 的 垂 向 运动 ,尽管 描述 了 土 体内 主要 的 水 
盐 运 动 过 程 ,但 缺乏 横向 和 各 向 异性 的 考虑 ,空间 
言 且 有 所 缺失 。 和 渗 过 程 中 的 优先 流 一 般 不 能 忽 
略 ,但 考虑 到 试验 区 地 面 平整 , 且 质 地 较 粗 砂 土 训 
面 均 一 ,没有 发 现 因 植物 根系 生长 或 者 动物 因素 而 
产生 的 孔隙 或 者 裂隙 ,发 生 优先 流 的 可 能 性 较 小 ， 
因而 本 次 试验 假定 优先 流 不 存在 。 此 外 在 模型 预 
测 时 ,没有 考虑 气候 变化 对 未 来 模型 输入 参数 的 影 
啊 , 后 续 研 究 工 作 会 着 重 考虑 此 方面 问题 。 


4 结论 

(1) 灌溉 后 ,表层 (0~30 cm) 土 壤 含水 率 迅 速 增 
加 ,可 快速 达到 1$%~20% 。 停 灌 后 又 快速 下 降 至 
10% 以 下 。50 em 以 下 的 土 层 含水 量变 化 较 小 ,通常 
在 5%~10% 之 间 波 动 。 
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要 好 于 下 层 土壤 
值 与 实测 值 的 空 
K , BURT te REREN 


(2) 对 于 土壤 水 分 的 数值 模拟 ,30 cm 以 上 深度 
;对 于 土壤 盐分 的 数值 模拟 ,模拟 
oe. 致 ,差异 性 不 

条 件 下 的 盐分 模拟 结果 可 信和 度 


更 高 。 


(3) 利用 HYDRUS-1D 模 型 对 高 矿 化 度 水 滴灌 


条 件 下 ,土壤 含水 率 和 含 盐 量 在 垂直 方向 的 变化 趋 


势 进行 


了 了 模拟 ,模拟 值 与 实测 值 拟 合 度 较 高 ,在 人 研 


究 区 的 应 用 性 较 强 ,能 够 为 沙漠 防护 林 的 建设 提供 
BAR SEF 
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Numerical simulation of water and salt migration in desert Soil in 


the lower reaches of Tarim River under salt-water irrigation 


WANG Shiming'***, FAN Jinlong'’**, ZHAO Ying’, ZHANG Taotao'***, LI Shengyu'’** 


(1. State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, 
Urumqi 830011, Xinjiang, China; 2. Mosuowan Desert Research Station, Xinjiang Institute of Ecology and Geography, Chinese 
Academy of Sciences, Shihezi 832000, Xinjiang, China; 3. Taklimakan Desert Research Station, Xinjiang Institute of Ecology and 
Geography Chinese Academy of Sciences, Korla 841000, Xinjiang, China; 4. University of Chinese Academy of Sciences, Beijing 
100049, China; 5. Ludong University, Yantai 264039, Shandong, China) 


Abstract: During shelterbelt construction in the lower reaches of Tarim River, Xinjiang, China, the water and salt 
balance of the soil profile clearly changed because of underground saltwater irrigation. Establishing a numerical 
model is a very efficient means of understanding the vertical water and salt dynamics of soil. By reconstructing 
the dynamic process of water and salt in the soil profile, the regularity of soil water and salt migration in a high- 
salt environment can be revealed, which would be useful for the construction of shelterbelts in ecological 
engineering projects and sustainable management of the subsequent irrigation. In this study, the distribution 
characteristics and temporal changes of the soil water and salt under different salinity conditions for groundwater 
irrigation in the early construction stages of the HYDRUS-1D shelter forest were simulated. Soil water and salt data 
were collected for optimization of the soil water characteristic parameters, and the applicability of the established 
numerical model was evaluated. The results showed that the water and salt contents of the surface soil (0-30 cm 
depth) were greatly affected by irrigation with large fluctuations. Additionally, the salt content was highly 
concentrated at the soil surface. In the deep layer (50-150 cm), irrigation had a small influence with small fluctuations 
in values. The simulated values of the numerical model were in good agreement with the measured values, which 
demonstrated that the model could reflect the soil water and salt migration at the test site. The established model 
can be used to simulate the soil water and salt migration in the desert area downstream of Tarim River under 
different salinity conditions for irrigation and provides a theoretical basis and technical support for the 
construction of an ecological shelterbelt in this area. Using the HYDRUS numerical simulation method 
effectively reduces the research costs of accurately predicting future changes in the soil water and salt 
environment, and it can also play an effective auxiliary role in research on other soil processes. This experiment 
was conducted in a desert environment, and few experimental studies have been performed in similar 
environments. Hence, this study is innovative and can provide new ideas and methods for experimental studies in 
relevant desert environments. 


Key words: salt water irrigation; soil water and salt migration, HYDRUS-1D model; lower Tarim River 


